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APRS  PERFORMANCE  EVALUATION  TESTS 


by 

R.  T.  Moore,  R.  Michael  McCabe  and  R.  Allen  Wilkinson 


INTRODUCTION 

The  FBI’s  Automatic  Fingerprint  Reader  Systems  (AFRS)  are  designed  to  scan  life-sized 
fingerprint  images  and  to  detect  selected  features  including,  but  not  always  limited  to, 
minutiae.  The  position  and  orientation  of  these  minutiae  are  recorded  in  units  of  X,  Y,  and 
Theta.  For  a variety  of  reasons,  an  AFRS  generally  fails  to  detect  all  of  the  true  minutiae 
in  a fingerprint  and  generally  makes  a number  of  false  detections.  Also,  the  true  minutiae 
that  are  detected  in  different  readings  of  the  same  fingerprint  impression  are  not 
necessarily  found  in  exactly  the  same  relative  positions.  Thus,  the  accuracy,  consistency 
and  reliability  of  minutiae  detection  are  useful  measures  of  reader  performance. 


PERFORMANCE  METRICS 

Historically,^  reader  performance  has  been  evaluated  by  superimposing  a plot  of  the 
minutiae  detections  made  by  the  fingerprint  reader  on  a photograph  of  the  fingerprint 
enlarged  to  the  same  scale  as  the  plot.  Then,  a fingerprint  expert  marks  the  plot 
indicating  each  true,  false,  or  missing  minutiae  in  the  area  of  the  fingerprint  that  is  being 
considered.  Limits  of  positional  accuracy  and  angular  accuracy  are  applied  and  each  minutia 
detected  by  the  reader  is  classed  as  either  a true  or  a false  minutia.  A detection  score  is 
then  assigned: 

Detection  Score  = (D  - M - F)  / T (1) 

where  D = number  of  true  minutiae  detected  by  reader 
M = number  of  true  minutiae  missed  by  reader 
F = number  of  false  minutiae  detected  by  reader 
T = total  number  of  true  minutiae  in  fingerprint 

It  should  be  noted  that  the  scoring  in  Equation  1 penalizes  a missed  minutia  twice  as 
heavily  as  a false  minutia. 

Scores  can  range  from  1.00  for  perfect  performance  (no  missed  or  false  minutiae)  down  to 
large  negative  values  which  might  result  with  few  true  detections  and  many  false  detections. 

This  procedure  for  evaluating  reader  performance  is  useful  but  quite  labor  intensive.  As  a 
consequence,  reader  performance  evaluation  in  the  past  has  usually  been  limited  to  the 


^ R.  M.  Stock  emd  C.  W.  Swonger,  Development  and  Evaluation  of  a Reader  of 
Fingerprint  Minutiae,  Cornell  Aeronautical  Laboratory  Technical  Report  CAL  No.  XM-2478-X- 
1,  Contract  J-FBI-6499,  January  1969. 
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examination  of  only  a small  number  of  sample  fingerprints.  Until  fairly  recently,  the 
measurement  of  reader  consistency,  or  repeatability,  has  received  very  little  attention. 

Recent  work  has  been  directed  toward  automating  certain  portions  of  the  process.  In 
particular,  computer  programs  have  been  developed^  that  permit  two  sets  of  minutiae  data  to 
be  aligned  in  translation  and  rotation  to  a "best  fit"  position.  These  programs  are  adapted 
from  the  M-40  matcher  algorithm  and  are  applied  in  four  stages  in  an  iterative  manner. 
From  the  final  positioning,  tolerances  are  applied  to  the  relative  displacement  and 
orientation  of  minutiae  pairs  that  are  nearly  aligned  with  each  other  to  determine  whether 
or  not  they  are  mates. 

If  one  of  the  two  sets  of  minutiae  data  is  considered  "ground  truth"  data,  then  the  number 
of  minutiae  pairs  that  are  within  the  tolerance  limits  establishes  the  value  of  D in  Equation 
1.  The  total  number  of  minutiae  in  the  "ground  truth"  data  set  establishes  the  value  of  T. 
The  value  of  M is  equal  to  T - D,  while  F is  equal  to  the  total  number  of  minutiae  detected 
in  the  data  set  being  compared  minus  the  value  of  D. 


COMPARISON  WITH  "GROUND  TRUTH" 

In  a typical  instance,  a fingerprint  is  read  manually  using  any  one  of  several  different 
types  of  semi-automatic  terminals.  The  reading  is  performed  twice;  once  by  each  of  two 
different  operators.  Next,  the  minutiae  from  each  reading  are  plotted  out  at  lOX 
enlargement.  The  two  plots  are  then  overlaid  and  any  discrepancies  are  easily  identified 
and  can  be  resolved.  The  resulting  minutiae  list,  after  all  anomalies  have  been  resolved,  is 
considered  a reasonable  estimate  of  "ground  truth"  for  that  fingerprint.  It  may  not 
exactly  coincide  with  the  machine  read  data  for  mating  pairs  of  minutiae  on  that 
fingerprint  simply  because  the  human  expert  and  the  machine  sometimes  use  slightly 
different  rules  in  assigning  position  and  orientation  to  a minutia,  but  usually  these 
differences  tend  to  be  small. 

When  the  data  from  a machine  read  fingerprint  is  compared  with  the  "ground  truth"  data, 
the  values  for  entry  in  Equation  1 can  be  developed.  In  addition,  the  distribution  of  the 
displacements  in  the  position  of  mating  minutiae  pairs  can  be  determined.  The 
displacement  is  calculated  as  the  straight  line  distance  between  each  of  the  mating  pairs 
whose  X,  Y and  Theta  differences  are  within  the  tolerance  limits  that  have  been 
established.  It  is  calculated  as: 

S =7(  x2  + y2)  (2) 

where  S = straight  line  distance 
X = displacement  in  X 
Y = displacement  in  Y 


When  "ground  truth"  minutiae  data  are  recorded,  each  minutia  is  given  a reference  number. 
Each  minutia  of  machine  read  minutiae  data  is  also  assigned  a reference  number.  The 
comparison  programs  use  these  reference  numbers  to  identify  mating  pairs  of  minutiae.  These 
reference  numbers  are  also  used  in  recording  how  frequently  each  "ground  truth"  minutiae  is 
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A description  and  listing  of  these  programs  appears  in  the  Appendix. 


2 


detected  on  successive  readings  of  the  same  fingerprint,  or  on  repeated  processing  of  the 
gray  scale  data  developed  during  a single  reading  of  the  fingerprint. 


COMPARISON  WITH  "FALSE  GROUND  TRUTH" 

The  programs  also  provide  an  alternative  means  for  evaluating  the  consistency  of  reader 
performance  given  identical  gray  scale  fingerprint  information  as  the  input  data  on 
successive  runs.  Here,  a "false  ground  truth"  data  set  is  calculated  and  used  to  determine 
how  repeatedly  each  minutia,  either  true  or  false,  is  detected. 

The  recorded  gray  scale  data  is  input  to  the  reader  the  first  time  and  the  minutiae  are 
detected.  Each  minutia  in  the  first  set  of  detections  from  the  fingerprint  is  considered  a 
potential  site  for  a "false  ground  truth"  minutia  "cluster"  location.  Then,  the  recorded 
gray  scale  data  is  input  a second  time  and  the  minutiae  are  detected  again.  The  minutiae  of 
the  second  set  of  detections  are  then  translated  and  rotated  for  "best  fit"  with  the  first  set 
of  minutiae.  The  position  and  orientation  of  the  minutiae  in  the  second  set  of  detections 
are  compared  with  those  of  the  first  and  those  that  are  within  the  selected  tolerance  limits 
are  declared  mates.  Each  mating  pair  of  minutiae  is  the  basis  for  adjustment  in  the 
location  of  its  "cluster"  site.  The  new  location  of  the  "cluster"  site  is  the  mean  position  of 
its  two  members.  Minutiae  that  are  beyond  the  tolerance  limits  for  displacement  and/or 
orientation  do  not  have  mates.  Their  locations  establish  candidate  new  sites  for  additional 
"clusters"  to  be  formed  with  minutiae  detected  from  the  third  and  subsequent  passes  of  the 
recorded  gray  scale  data  through  the  reader. 

This  process  is  repeated  on  the  subsequent  minutiae  detections  made  from  the  fingerprint 
data.  As  additional  minutiae  on  the  subsequent  runs  are  found  to  be  located  within  the 
tolerance  limits  of  position  and  orientation  that  are  established  for  mates,  they  become 
members  of  that  "cluster",  or  become  potential  sites  for  new  "clusters".  The  center  point 
of  each  "cluster"  is  continually  recalculated  on  the  basis  of  the  positions  of  the  members  of 
that  "cluster".  Occasionally,  minutiae  will  leave  one  "cluster"  and  join  another  as  a result 
of  this  recalculation  of  center  position  of  the  "cluster". 

An  analysis  of  these  "cluster"  sizes  provides  an  indication  of  reader  consistency.  Ideally, 
all  "clusters"  should  have  a number  of  members  equal  to  the  number  of  times  that  the 
identical  gray  scale  data  was  passed  through  the  system  for  minutiae  detection.  Smaller 
"clusters"  are  indicative  of  inconsistencies  in  performance. 


COMPARISON  WITH  FIRST  PASS  DATA 

A third  means  of  evaluating  consistency  is  provided  when  minutiae  data  from  the  first  of  a 
series  of  passes  of  gray  scale  data  are  used  as  a reference.  These  data  are  compared  with 
the  minutiae  data  from  each  of  the  succeeding  passes  of  the  identical  gray  scale  data.  The 
"cluster"  size  information  developed  from  this  routine  is  similar  but  not  identical  to  that 
developed  from  use  of  the  "false  ground  truth"  procedure  described  above.  This  is  because 
there  is  no  recalculation  of  "cluster"  center  position  as  members  are  added  to  it.  In 
addition  to  cluster  size  information,  this  routine  also  provides  information  on  the 
distribution  of  the  values  of  the  displacement  in  the  position  of  the  mating  minutiae.  Since 
the  same  gray  scale  information  is  being  input  on  each  pass,  this  displacement  is  assumed  to 
be  caused  by  noise  in  the  minutiae  detection  electronics. 


3 


TEST  RESULTS 


In  connection  with  the  conversion  of  the  APRS  from  a flying  spot  scanner  to  a solid  state 
scanner,  certain  special  test  materials  were  prepared.  Among  these  was  a recording  of  the 
fingerprints  of  a single  individual  using  three  different  inking  densities,  "Light",  "Medium", 
and  "Heavy".  These  fingerprints  were  scanned  by  the  digital  scainner  on  APRS  No.  4 and 
the  gray  scale  data  were  recorded  on  disk.  This  permitted  the  same  gray  scale  data  to  be 
directed  repeatedly  to  the  APRS  preprocessor  for  minutiae  detection.  Ideally,  the  same  gray 
scale  data  should  produce  cm  identical  minutiae  list  each  time  it  is  processed.  Differences 
in  the  identity  or  position  of  the  detected  minutiae  represent  imperfect  performance  of  the 
preprocessor  which  might  be  caused  by  noise  or  other  factors. 

This  gray  scale  data  recorded  on  disk  was  input  to  the  preprocessor  of  APRS  No.  4 eight 
times.  The  detected  minutiae  were  registered  and  clipped  and  then  recorded  as  eight 
fingers  of  data  from  each  of  three  cards.  The  data  from  the  lightly  inked  card  were 
recorded  as  card  105,  medium  as  card  203,  and  heavy  as  card  301  from  system  No.  4. 

The  same  procedure  was  followed  using  the  same  fingerprints  on  APRS  No.  2.  This  system 
has  had  decoupling  capacitors  installed  on  the  printed  circuit  cards  in  the  preprocessor. 


"GROUND  TRUTH-  MATCH  RESULTS 

The  "ground  truth"  data  for  these  three  fingerprints  was  obtained  using  the  Graphic  Pen^ 
at  NBS.  The  objective  was  to  record  only  those  minutiae  that  appeared  within  the  area 
covered  by  the  clipping  box.  This  would  provide  minutiae  from  an  area  that  was  common  to 
the  area  of  the  machine-read  minutiae  data.  A plot  of  the  machine  read  data  was  centered 
in  the  field  of  view  of  the  Graphic  Pen  in  order  to  approximately  define  the  clipping  box 
boundaries.  Then  the  fingerprint  was  positioned  under  the  magnifier  lenses  and  adjusted  to 
make  the  minutiae  agree  with  those  on  the  machine-read  plot.  Next  the  outline  of  the 
clipping  box  was  centered  on  the  machine-read  plot.  Finally  the  minutiae  of  the  fingerprint 
that  were  within  or  even  slightly  beyond  the  boundaries  of  the  clipping  box  were  read  by 
two  different  operators  to  develop  the  "ground  truth"  data  base  for  the  minutiae  of  that 
finger. 

There  may  be  an  element  of  uncertainty  in  this  process.  This  comes  about  because  in  the 
machine-read  data,  the  clipping  box  is  centered  on  a position  that  is  found  as  a result  of 
processing  the  ridge  flow  data  derived  from  the  machine-read  gray  scale  information.  If 
noise  in  the  system  affects  minutiae  positions,  it  may  also  affect  ridge  flow  data.  This 
might  cause  the  clipping  box  boundaries  to  be  slightly  different  for  each  pass  of  identical 
gray  scale  data  through  the  preprocessor.  Because  of  this,  some  of  the  minutiae  near  the 
boundaries  of  the  clipping  box  in  the  "ground  truth"  data  set  may  actually  be  outside  the 
boundaries  of  the  clipping  box  for  the  machine-read  data  on  one  or  more  of  the  processing 
passes.  Candidate  minutiae  have  been  identified  that  might  have  been  missed  for  this 
reason.  These  candidates  consist  of  "ground  truth"  minutiae  that  were  located  very  close 


R.  T.  Moore  and  J.  R.  Park,  "The  Graphic  Pen  - An  Economical  Semiautomatic 
Fingerprint  Reader",  Proc.  1977  Carnahan  Conference  on  Crime  Countermeasures,  UKY 
BU112,  ORES,  Univ  of  Kentucky,  Lexington,  KY  40506. 
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to  the  clipping  box  boundary  that  were  not  detected  in  any  of  the  multiple  passes  of  the 
gray  scale  data  through  the  preprocessor. 

On  the  lightly  inked  fingerprint  there  were  83  minutiae  on  the  initial  "ground  truth"  list. 
Of  these,  nine  minutiae,  Nos.  1,  2,  3,  4,  8,  12,  23,  56,  and  78  were  not  detected  in  any 
pass  of  the  data.  It  is  reasonable  to  assume  that  they  might  have  been  outside  the 
clipping  box.  A value  of  74  is  therefore  used  for  T in  Equation  1 for  this  fingerprint.  On 
the  fingerprint  with  medium  inking,  there  were  75  minutiae  on  the  initial  list  and  three  of 
these.  Nos.  1,  49  and  53  were  not  detected  in  any  pass  of  the  data,  so  T is  assumed  to  have 
a value  of  72.  On  the  fingerprint  with  the  heavy  inking,  there  were  76  minutiae  on  the 
initial  "ground  truth"  list.  Four  of  these.  Nos.  5,  6,  19  and  51  had  no  detections  in  any 
pass  of  the  machine-read  data  so  a value  of  72  is  assumed  for  T. 

The  tolerance  limits  on  displacement  and  orientation  for  a pair  of  minutiae  to  be 
considered  mates  are  0.4  mm  (four  X,  Y matcher  units)  and  12  degrees  (Theta). 

With  the  gray  scale  data  from  these  three  fingerprints  recorded  on  disk  and  then  entered 
as  eight  independent  passes  into  the  preprocessor  of  APRS  No.4,  the  values  for  use  in 
Equation  1 are  shown  in  the  following  Table  1. 


TABLE  1 

Detection  Scores  for  APIS  No.  4 


Lieht 

Medium 

Heaw 

Pass 

D 

M 

F 

D 

M 

F 

D 

M 

F 

1 

18 

56 

60 

20 

52 

54 

28 

44 

60 

2 

20 

54 

46 

28 

44 

49 

30 

42 

57 

3 

14 

60 

48 

27 

45 

43 

35 

37 

58 

4 

20 

54 

51 

19 

51 

47 

34 

38 

49 

5 

18 

56 

57 

24 

48 

51 

33 

39 

53 

6 

15 

59 

57 

22 

50 

40 

40 

32 

46 

7 

16 

58 

51 

24 

48 

53 

33 

39 

54 

8 

19 

55 

50 

24 

48 

45 

34 

38 

49 

Tot. 

452 

420 

188 

388 

382 

267  309 

426 

Score;  -732/592  = -1.24  -582/576  = -1.01  -468/576  = -0.81 

The  mean  number  of  true  minutiae  detected  on  the  lightly  inked  print  was  17.5  or  24%  of 
the  "ground  truth"  minutiae.  On  the  medium  inked  print  it  was  23.5,  or  33%  of  "ground 
truth".  On  the  heavily  inked  print  it  was  33.375,  or  46%  of  "ground  truth". 

On  the  lightly  inked  print,  34  different  true  minutiae  were  detected  on  one  or  more  passes, 
but  only  six  were  consistently  detected  on  all  eight  passes.  On  the  medium  inked  print,  46 
different  true  minutiae  were  detected  on  one  or  more  passes,  and  only  five  were  detected 
on  all  eight  passes.  With  the  heavily  inked  print,  48  true  minutiae  were  detected  on  one  or 
more  passes  and  16  were  detected  on  all  eight  passes.  The  distribution  of  the  numbers  of 
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minutiae  and  consistency  of  their  detection  with  eight  passes  through  the  preprocessor  of 
APRS  No.  4 is  shown  in  Table  2. 


TABLE  2 


Consistency  of  Detections.  APRS  No,  4 
LIGHT  MEDIUM  HEAVY 


7 

5 

6 

3 
2 
1 

4 
6 


11 

5 

8 

2 

3 

7 

5 

5 


3 

3 

5 

6 
6 
1 
8 

16 


Minutiae  detected  1 time 
•’  " 2 times 


3 

4 

5 

6 

7 

8 


34 


46  48  Minutiae  detected  one  or  more  times 


Pigure  1 shows  this  same  data  in  a different  way.  This  is  a cumulative  distribution  of  the 
percentage  of  the  minutiae  detections  that  are  repeated  eight  times,  seven  times,  etc. 

Pigures  2,  3 and  4 show  distribution  of  distances  that  the  minutiae  were  displaced  from  the 
"ground  truth"  position  in  each  of  the  eight  passes  of  the  gray  scale  data  through  the  APRS 
No.  4 preprocessor  for  the  light,  medium  and  heavily  inked  image.  Pigure  5 shows  a 
summary  of  this  same  information  for  the  eight  passes  and  all  three  inkings.  This  figure 
shows  that  the  most  frequent  displacement  is  two  units  with  the  lightly  inked  images.  The 
heavily  inked  image  shows  a noticeable  peak  at  one  unit  of  displacement. 

The  preprocessor  of  APRS  No  2.  has  been  modified  with  the  addition  of  some  decoupling 
capacitors  to  its  printed  circuit  cards  in  an  attempt  to  reduce  internally  generated  noise. 
The  same  procedures  were  followed  with  this  reader  as  with  APRS  No.  4.  Gray  scale  data 
was  recorded  on  disk  and  passed  through  the  preprocessor  eight  times  and  the  minutiae 
detections  from  each  pass  compared  with  the  "ground  truth"  data  from  the  light,  medium  and 
heavily  inked  fingerprint.  The  results  are  shown  in  Table  3.  They  are  very  comparable  to 
the  performance  shown  for  APRS  No.  4 in  Table  1.  The  lightly  inked  print  provided  a mean 
of  17.5  true  minutiae  or  24%  of  "ground  truth".  On  the  medium  inked  print  the  mean 
number  of  detections  was  21.625  or  30%  of  "ground  truth".  The  heavily  inked  print  yielded 
a mean  of  33.125  minutiae  which  is  46%  of  the  "ground  truth". 

The  distribution  of  numbers  of  minutiae  and  the  consistency  of  detection  is  shown  in  Table 
4. 
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TABLES 


Detection  Scores  for  APRS  No.  2 


LIGHT 

MEDIUM 

HEAVY 

Pass 

D 

M 

F 

D 

M 

F 

D 

M 

F 

1 

24 

50 

53 

18 

54 

55 

34 

38 

50 

2 

17 

57 

41 

28 

44 

57 

32 

40 

56 

3 

20 

54 

51 

21 

51 

53 

31 

41 

60 

4 

16 

58 

47 

23 

49 

52 

36 

36 

45 

5 

12 

62 

48 

20 

52 

41 

37 

35 

46 

6 

21 

53 

62 

21 

51 

57 

31 

41 

54 

7 

15 

59 

58 

21 

51 

51 

31 

41 

54 

8 

15 

59 

62 

21 

51 

55 

32 

40 

58 

Tot. 

452  422 

m 

403 

4n 

2^ 

430 

Score: 

-733/592  = 

-1.24 

-64S/576  = 

-1.13 

-476/576 

= -0.83 

TABLE  4 

Consistencv  of  Detections.  AFRS  No.  2 

LIGHT 

MEDIUM 

HEAVY 

4 

11 

6 Minutiae  detected  1 time 

8 

4 

2 

" 2 times 

2 

3 

3 

" 3 

1 

11 

4 

" 4 '• 

2 

2 

2 

" 5 " 

5 

4 

5 

" 6 " 

6 

4 

7 

If  ^ II 

3 

5 

17 

" 8 " 

W 

44” 

46  Minutiae  detected  one  or  more  times 

Figure  6 shows  the  cumulative  distribution  of  the  repeated  minutiae  detections  shown  in 
Table  4. 

Figures  7,  8 and  9 show  distribution  of  distances  that  the  minutiae  were  displaced  from  the 
"ground  truth"  position  in  each  of  the  eight  passes  of  the  gray  scale  data  through  the  AFRS 
No.  2 preprocessor  for  the  light,  medium  and  heavily  inked  image. 

Figure  10  shows  a summary  of  this  same  information  for  the  eight  passes  and  all  three 
inkings.  This  figure  shows  that  the  most  frequent  displacement  is  two  units  with  the  light 
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and  medium  inked  images.  The  heavily  inked  image  again  shows  a noticeable  peak  at  only 
one  unit  of  displacement. 


"FALSE  GROUND  TRUTH"  RESULTS 

"False  ground  truth"  is  the  name  that  has  been  given  to  the  set  of  candidate  "cluster"  sites 
in  (X,  Y,  Theta)  space  representing  every  minutiae  detection,  true  or  false,  resulting  from 
multiple  passes  of  the  same  gray  scale  data  through  the  preprocessor  and  minutiae 
detection  logic.  Each  field  of  detected  minutiae  is  translated  and  rotated  to  a "best  fit" 
position  with  respect  to  the  field  of  minutiae  detected  in  the  first  pass.  Mating  minutiae 
are  identified  and  the  "cluster"  site  is  recalculated  to  be  the  mean  position  and  orientation 
of  the  mates  comprising  the  cluster.  Minutiae  that  do  not  have  mates  within  the  tolerance 
limits  of  displacement  and  rotation  still  establish  candidate  sites  that  may  become  populated 
with  subsequent  passes  of  the  data. 

Since  the  gray  scale  data  that  is  used  is  digital,  performance  in  a noise-free  environment 
would  be  expected  to  produce  N "clusters"  each  having  a population  of  M mates,  where  N 
is  a constant  number  of  minutiae  detections  and  M is  the  number  of  times  that  the  same 
gray  scale  data  is  passed  through  the  detector.  The  fact  that  this  is  not  the  case  is  a 
matter  of  serious  concern. 

Table  5 shows  the  performance  of  AFRS  No.  4 on  the  light,  medium  and  heavily  inked 
prints.  The  same  data  is  shown  in  Fig.  11  as  the  cumulative  distribution  of  the  percentage 
of  multiple  detections.  There  is  very  small  difference  in  the  results  from  the  light  and 
medium  inked  print,  while  the  heavily  inked  print  produced  noticeably  more  consistent 
detections. 

Table  6 and  Fig.  12  show  the  corresponding  performance  of  AFRS  No.  2.  Here  the 
superiority  of  the  heavily  inked  print  is  much  less  pronounced. 


TABLE  5 

SUMMARY  OF  "FALSE  GROUND  TRUTH"  CLUSTER  SIZES  - AFRS  NO.  4 


LIGHT  MEDIUM  HEAVY 


60 

69 

33 

Minutiae  detected  1 time 

27 

30 

16 

" " 2 times 

21 

16 

9 

" " 3 

13 

12 

19 

" 4 " 

16 

8 

11 

" " 5 " 

6 

14 

7 

" 6 " 

9 

11 

12 

tl  H ^ M 

19 

18 

43 

" " 8 " 

itT 

150 

Different  minutiae  detected 

one  or  more  times 
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TABLE  6 


SUMMARY  OF  "FALSE  GROUND  TRUTH"  CLUSTER  SIZES  - APRS  N0.2 

LIGHT  MEDIUM  HEAVY 


58 

53 

57 

Minutiae  detected  1 time 

36 

27 

17 

" " 2 times 

14 

11 

12 

" 3 

15 

14 

13 

" " 4 " 

8 

9 

8 

" " 5 " 

8 

8 

16 

" 6 " 

14 

15 

12 

" " 7 " 

18 

25 

37 

" " 8 " 

m 

Different  minutiae  detected 

one  or  more  times 


The  software  routines  that  developed  the  data  shown  above  in  Tables  5 and  6 for  "false 
ground  truth"  performance  hst  the  identity  of  each  of  the  minutiae  that  formed  each  of 
the  clusters.  Many  of  these  clusters  appear  to  have  been  formed  from  false  detections 
that  occurred  repeatedly.  A comparison  of  the  numbers  of  minutiae  detected  eight  times  in 
Tables  2 and  4 with  the  number  of  minutiae  detected  eight  times  in  Tables  5 and  6 reveals 
information  about  how  frequently  these  repeated  false  detections  occurred. 


FIRST  PASS  MATCH  RESULTS 

This  comparison  takes  the  minutiae  detections  resulting  from  the  first  pass  of  the  gray  scale 
data  from  the  light,  the  medium,  and  the  heavily  inked  fingerprint  and  treats  it  as  the 
reference  against  which  the  other  seven  passes  of  the  data  from  each  of  these  fingerprints 
are  compared.  As  is  the  case  of  the  "false  ground  truth"  data  there  is  no  distinction  made 
as  to  whether  the  minutiae  detected  are  true  or  false.  The  chief  utility  of  this  comparison 
is  to  provide  a measure  of  the  displacements  in  position  and  orientation  of  mating  minutiae 
that  is  not  biased  by  the  human  position  selection  rules  as  it  is  in  the  "ground  truth"  data. 

Figures  13  through  15  show  the  seven  individual  results  obtained  from  each  of  the  three 
degrees  of  inking  on  APRS  No.  4.  Figures  16  and  17  show  the  summary  data  for  this 
system. 

Figures  18  through  20  show  the  seven  individual  results  from  each  of  the  three  inkings  on 
AFRS  No.2  and  Figs.  21  and  22  cover  the  summary  results  for  this  system. 

These  data  show  that  the  most  frequently  observed  displacement  in  position  for  mating 
minutiae  is  one  matcher  unit,  although  there  are  displacements  greater  than  that  in  a few 
percent  of  the  detections.  With  no  noise,  it  would  be  expected  that  there  would  be  no 
displacement,  since  the  same  digital  data  is  input  each  time. 
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CONCLUSIONS 


These  test  results  indicate  that  despite  the  improvements  resulting  from  upgrading  the 
readers  with  the  new  sohd  state  scanner  subsystems,  there  are  serious  problems  in  the 
APRS.  These  are  manifest  in  the  form  of  inconsistencies  in  the  detection  of  minutiae, 
both  true  and  false,  even  when  identical  gray  scale  data  is  input  to  the  preprocessor  and 
minutiae  detection  circuitry.  The  performance  is  more  erratic  when  lightly  or  moderately 
inked  fingerprints  are  used  than  with  heavily  inked,  high  contrast,  images.  In  these  tests, 
only  about  one  quarter  of  the  true  minutiae  were  detected  in  lightly  inked  prints,  one  third 
in  the  moderately  inked  prints  and  one  half  in  the  heavily  inked  prints  on  any  single 
processing  of  the  fingerprint  data.  Presumably  the  observed  inconsistencies  in  performance 
result  from  noise  in  the  preprocessor.  Minor  inconsistencies  in  detection  performance  could 
be  expected  to  be  caused  by  the  recursive  behavior  of  the  ridge  valley  filter,  but  these 
would  be  expected  to  be  constrained  to  the  top  few  per  cent  of  the  fingerprint  image  area. 
Since  several  of  the  instances  where  minutiae  were  detected  consistently  (eight  times)  occur 
in  the  top  five  percent  of  the  image  area,  it  is  not  believed  that  the  recursive  attributes  of 
the  filter  are  a major  contributor  to  the  inconsistencies  observed  in  minutiae  detection. 

It  is  believed  that  the  detection  probability  displayed  in  these  tests  is  not  capable  of 
supporting  an  effective  automated  latent  fingerprint  identification  system.  It  is  suggested 
that  long  range  planning  should  contemplate  reconversion  of  the  files  to  be  used  for  an 
automated  latent  system.  Rigorous  quality  control  measures  are  suggested  to  insure  that 
the  quality  of  the  re-converted  file  data  is  maintained  at  an  acceptable  level.  It  is 
believed  that  some  of  the  software  routines  listed  in  the  Appendix  would  be  appropriate 
candidates  for  use  in  support  of  this  function. 

These  test  data  also  strongly  suggest  that  there  is  no  significant  difference  in  the 
performance  of  APRS  No.  2 and  No.  4 and  that  the  FBI  made  the  correct  choice  in  not 
adding  the  decoupling  capacitors  to  the  remainder  of  the  APRS  preprocessors.  It  is 
believed  that  much  more  extensive  measures  will  be  required  to  correct  the  problems  in 
these  systems. 

Finally,  these  data  show  a definite  trend  to  improved  performance  with  increasing  ink 
density.  This  is  in  agreement  with  vendor  claims  that  readers  perform  better  with 
electronically  generated  fingerprints  which  provide  images  with  high  contrast. 
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CUM.  % BY  NO.  OF  DETECTIONS 


Figure  1 
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CUMULATIVE  % OF  TOTAL  DETECTIONS  - GROUND  TRUTH 
SAME  GRAY  SCALE  DATA  - AFRS  |4  - ALL  3 INKINGS 


NO  OF  MINUriAC  MO  OF  UINUTIAE  NO  OF  UINUTIAE  NO.  OF  MINUTIAE 


MINUTIAE  DISPLACEMENT  FROM  GROUND  TRUTH 
SAME  GRAY  SCALE  DATA  - APRS  #4  - LIGHT  INKING 


Figure  2 
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NO.  OF  MINUTIAE  NO  OF  MINUTIAE  NO.  OF  MINUTIAE  HO  OF  MINUTIAE 


MINUTIAE  DISPLACEMENT  FROM  GROUND  TRUTH 
SAME  GRAY  SCALE  DATA  - APRS  #4  - MEDIUM  INKING 


Figure  3 
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NO  OF  MINUr(AE  NO  OF  MINUTIAE  NO.  OF  MINUTIAE  NO.  OF  MINUTIAE 


MINUTIAE  DISPLACEMENT  FROM  GROUND  TRUTH 
SAME  GRAY  SCALE  DATA  - AFRS  #4  - HEAVY  INKING 


Figure  4 
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STAR;  LIGHT  TRIANGLE;  MEDIUM  CIRCLE;  HEAVY 


NO.  OF  MINUTIAE 


Figure  5 
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MINUTIAE  DISPLACEMENT  FROM  GROUND  TRUTH 
SAME  GRAY  SCALE  DATA  - AFRS  |4  - ALL  3 INKINGS 


CUM.  % BY  NO.  OF  DETECTIONS 


Figure  6 
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CUMULATIVE  % OF  TOTAL  DETECTIONS  - GROUND  TRUTH 
SAME  GRAY  SCALE  DATA  - AFRS  |2  - ALL  3 INKINGS 


NO  OF  MINUTIAE  NO  OF  MINUTIAE  NO.  OF  MINUTIAE  NO.  OF  MINUTIAE 


MINUTIAE  DISPLACEMENT  FROM  GROUND  TRUTH 
SAME  GRAY  SCALE  DATA  - AFRS  #2  - LIGHT  INKING 


figure  7 
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NO.  OF  MINUTIAE  NO  OF  MINUTIAE  NO.  OF  MINUTIAE  NO  OF  MINUTIAE 


MINUTIAE  DISPLACEMENT  FROM  GROUND  TRUTH 
SAME  GRAY  SCALE  DATA  - APRS  #2  - MEDIUM  INKING 


Figure  8 
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NO  OF  UINUriAE  NO  OF  UINUTIAE  NO.  OF  MINUIIAE  NO.  OF  MINUTIAE 


MINUTIAE  DISPLACEMENT  FROM  GROUND  TRUTH 
SAME  GRAY  SCALE  DATA  - APRS  §2  - HEAVY  INKING 


Figure 


19 


9 


STAR;  LIGHT  TRIANGLE:  MEDIUM  CIRGLE:  HEAVY 


NO.  OF  MINUTIAE 


Figure  10 
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MINUTIAE  DISPLACEMENT  FROM  GROUND  TRUTH 
SAME  GRAY  SCALE  DATA  - AFRS  §2  - ALL  3 INKINGS 


CUM.  % BY  NO.  OF  DETECTIONS 


Figure  11 
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CUMULATIVE  % OF  TOTAL  DETECTIONS  - FALSE  GROUND  TRUTH 
SAME  GRAY  SCALE  DATA  - AFRS  |4  - ALL  3 INKINGS 


CUM.  % BY  NO.  OF  DETECTIONS 


Figure  12 
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CUMULATIVE  % OF  TOTAL  DETECTIONS  - FALSE  GROUND  TRUTH 
SAME  GRAY  SCALE  DATA  - AFRS  §2  - ALL  3 INKINGS 


MINUTIAE  DISPLACEMENT  FROM  PASS  #1 
SAME  GRAY  SCALE  DATA  - APRS  #4  - LIGHT  INKING 


Figure  13 
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NO  Of  MINUFIAE  NO.  OF  MINUTIAE  NO.  OF  MINUTIAE  NO  OF  MINUTIAE 


MINUTIAE  DISPLACEMENT  FROM  PASS  #1 
SAME  GRAY  SCALE  DATA  - AFRS  #4  - MEDIUM  INKING 


Figure  14 
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NO.  OF  MINUriAE  NO  OF  UINUIIAC  NO.  OF  UINUTIAC  NO.  OF  MINUTIAE 


MINUTIAE  DISPLACEMENT  FROM  PASS  #1 
SAME  GRAY  SCALE  DATA  - APRS  #4  - HEAVY  INKING 


Figure  15 
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STAR:  LIGHT  TRIANGLE:  MEDIUM  CIRCLE:  HEAVY 
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Figure  16 
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MINUTIAE  DISPLACEMENT  EROM  PASS  |1 
SAME  GRAY  SCALE  DATA  - APRS  |4  - ALL  3 INKINGS 


CUM.  % BY  NO.  OF  DETECTIONS 


Figure  17 
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CUMULATIVE  % OF  TOTAL  DETECTIONS  - WITH  PASS  1 
SAME  GRAY  SCALE  DATA  - APRS  |4  - ALL  3 INKINGS 


MINUTIAE  DISPLACEMENT  FROM  PASS  #1 
SAME  GRAY  SCALE  DATA  - APRS  #2  - LIGHT  INKING 


Figure  18 
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MINUTIAE  DISPLACEMENT  FROM  PASS  #1 
SAME  GRAY  SCALE  DATA  - AFRS  #2  - MEDIUM  INKING 


Figure 
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MINUTIAE  DISPLACEMENT  FROM  PASS  #1 
SAME  GRAY  SCALE  DATA  - APRS  #2  - HEAVY  INKING 


Figure 
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STAR:  LIGHT  TRIANGLE;  MEDIUM  CIRGLE:  HEAVY 
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Figure  21 
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MINUTIAE  DISPLACEMENT  FROM  PASS  |1 
SAME  GRAY  SCALE  DATA  - APRS  |2  - ALL  3 INKINGS 


CUM.  % BY  NO.  OF  DETECTIONS 


Figure  22 


32 


CUMULATIVE  % OF  TOTAL  DETECTIONS  - WITH  PASS  1 
SAME  GRAY  SCALE  DATA  - AFRS  §2  - ALL  3 INKINGS 


APPENDIX 


This  appendix  contains  a verbal  description,  flow  diagrams,  and  the  entire  code  for  the 
program  FINDTRANS  and  its  subroutines.  The  program  and  subroutines  are  written  in 
FORTRAN-77.  Only  two  of  the  subroutines,  FBIOPEN  and  FBIREAD  are  machine  specific  and 
only  run  on  a VAX/VMS  system,  while  the  remaining  subroutines  and  program  should  be 
portable. 

The  program  FINDTRANS  calculates  the  "best  fit"  transformation  values  for  delta  X,  delta  Y 
and  delta  THETA  to  be  used  in  matching  a pair  of  fingerprints.  Additionally,  the  program 
evaluates  the  accuracy  of  a match  performed  using  these  transformation  values.  A "best  fit" 
may  be  defined  as  the  orientation  in  translation  and  rotation  of  a comparison  fingerprint 
with  a mating  base  fingerprint,  such  that  the  positions  of  the  majority  of  minutiae  from  the 
comparison  fingerprint  are  relatively  close  to  those  of  the  base  fingerprint  as  evidenced  by 
highest  matching  score.  A "best  fit"  transformation  may  be  one  of  several  possible  "best  fit" 
transformations  because  of  the  many  combinations  of  transformations  that  can  be  generated 
with  equal  scores.  In  this  situation  the  first  "best  fit"  transformation  found  is  the 
transformation  used. 

As  input,  FINDTRANS  requires  the  name  of  the  file  containing  the  base  fingerprint  (the  one 
which  will  not  be  transformed)  and  the  quantity  of  other  prints  to  be  compared  against  this 
base  fingerprint.  For  each  comparison  print,  the  minutiae  positions  are  read  in  and  the 
initial  "best  fit"  is  assumed  at  the  transformation  values  of  0 in  X,  0 in  Y,  0 in  rotation. 
FINDTRANS  then  issues  four  calls  to  the  subroutine  BESTFIT,  which  calculates  new 
transformation  values.  Each  time  BESTFIT  is  called,  a constant  theta  difference  limit  of  12 
degrees  is  used.  However,  the  straight  line  distance  tolerance  decreases  with  each  call  to 
BESTFIT.  The  values  that  are  used  for  the  four  calls  are  30,  15,  8,  and  4 respectively.  Only 
four  calls  of  BESTFIT  are  used  as  the  "best  fit"  values  after  the  fourth  iteration  show  very 
little  sign  of  change.  When  FINDTRANS  is  done,  the  "best  fit"  transformation  is  known  and 
diagnostic  output  of  the  match  performed  using  this  transformation  is  generated. 

BESTFIT  is  the  primary  subroutine  called  by  FINDTRANS,  as  it  performs  the  majority  of  the 
processing.  It  initially  saves  the  assumed  "best  fit"  values  that  were  passed  to  it  and  asks 
the  user  to  input  the  beginning  and  ending  rotational  range  with  values  stated  in  degrees. 
This  range  is  used  to  test  for  the  best  transformation  rotation  value.  This  will  be  a value 
between  the  assumed  "best  fit"  rotation  plus  the  beginning  value,  increasing  in  increments  of 
one  degree,  to  the  assumed  "best  fit"  rotation  plus  the  ending  value.  The  beginning  and 
ending  values  may  be  negative  but  the  beginning  value  must  always  be  less  than  the  ending 
value.  The  comparison  print  is  transformed  using  the  present  rotation  value  of  the 
programming  loop.  After  the  transformation,  the  two  prints  are  compared  to  find  possible 
matches  in  minutiae  positions.  Any  minutiae  position  in  the  base  print  that  has  a minutiae 
position  in  the  comparison  print  located  within  the  straight  line  distance  tolerance  and  with 
the  difference  in  theta  values  less  than  or  equal  to  12  degrees  is  considered  to  be  a 
possible  match.  These  possible  matches  are  saved  in  a table  called  the  S-table.  Each  entry  in 
the  S-table  is  given  a score  which  rates  the  closeness  of  the  match  between  the  two 
minutiae.  Using  a scoring  strategy,  a score  is  generated  for  the  entire  transformation.  If 
this  is  the  first  transformation  then  the  score  is  assumed  to  be  the  highest  score  and  is 
saved  with  the  transformation  values  that  generated  the  score.  If  this  is  not  the  first 
transformation  then  the  score  is  compared  against  the  saved  score.  If  the  new  score  is 
higher  than  the  saved  score,  then  the  values  and  score  are  saved  as  the  new  "best  fit"  and 
high  score.  When  the  last  rotation  value  of  the  loop  has  been  processed,  the  highest  score 
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and  the  corresponding  transformation  values  are  saved.  This  rotation  value  is  used  to 
retransform  the  comparison  print  and  then  recreate  the  S-table  for  this  transformation.  The 
data  in  the  S-table  is  used  to  determine  where  the  concentration  of  S-table  entries  are 
located.  This  point  of  concentration  is  the  negative  offset  needed  for  the  "best  fit"  X and  Y 
values.  If  the  point  of  concentration  is  at  position  (3,-1)  then  the  new  "best  fit"  values  are 
the  previous  "best  fit"  X value  minus  3,  and  the  previous  "best  fit"  Y value  plus  1.  BESTFIT 
returns  to  FINDTRANS,  its  caller,  the  "best  fit"  values  it  has  generated  and  the  score  of 
the  match  generated  with  these  transformation  values. 

Several  scoring  strategies  have  been  tested  in  an  effort  to  evaluate  the  accuracy  of  a match 
between  a base  fingerprint  and  a transformed  comparison  fingerprint.  The  strategies  that 
were  tested  include  highest  S-table  entry  score,  average  of  the  S-table  entry’s  scores  and 
total  of  the  S-table  entry’s  scores.  The  highest  S-table  score  was  found  to  be  an  effective 
method  until  the  prints  had  at  least  one  minutiae  position  match  exactly  on  both  prints. 
This  caused  the  highest  possible  score  to  be  generated.  Many  transformation  values  caused 
this  situation.  Due  to  this,  the  highest  score  method  could  not  accurately  choose  a "best 
fit".  To  solve  this  problem  more  than  one  S-table  entry  value  had  to  be  used  so  the  total 
and  average  methods  were  tested.  Averaging  the  S-table  on  some  matches  cause  a selection 
of  two  or  more  "best  fit"  transformations.  With  examination,  these  were  not  found  to  be 
very  accurate.  By  using  the  total  method  the  best  results  were  obtained.  The  total  method 
summed  the  S-table  scores  to  generate  the  overall  score.  This  method  was  used  in  BESTFIT 
to  create  a match  score. 

Example  of  results  using  all  three  methods 


Entry  # 

S-table  1 scores 

S-table  2 scores 

S-table  3 s 

1 

100 

120 

180 

2 

100 

110 

100 

3 

50 

70 

80 

4 

60 

80 

40 

5 

150 

90 

60 

6 

20 

No  entry 

No  entry 

Total  method  480 

470 

460 

Average  method  80 

94 

92 

Highest  method  150 

120 

180 

To  create  the  S-table  for  a comparison,  the  routine  MAK40  is  called.  This  routine  creates  a 
table  of  pairs  of  minutiae  whose  delta  X,  delta  Y and  delta  theta  differences  are  within  the 
tolerance  levels.  The  tolerance  limits  are  linear  distance,  usually  no  more  than  four  units, 
and  angular  difference,  usually  no  more  than  12  degrees.  This  table  keeps  track  of  the 
minutiae  identity  numbers  for  the  transformed  print  and  base  print,  the  differences  in  X, 
and  Y,  the  delta  theta  of  these  two  minutiae,  and  the  score  for  the  entry.  When  an  entry  is 
added  to  the  list  it  is  assigned  a score  of  zero.  Once  all  the  entries  have  been  created, 
MAK40  generates  a score  for  each  one.  To  do  this,  the  linear  distance  between  these  two 
entries  is  calculated  using  their  X and  Y difference  values  [1].  If  the  distance  is  not  within 
the  straight  line  distance  tolerance  level  or  is  greater  than  five,  it  has  no  effect  on  the 
score.  If  the  distance  is  within  the  straight  line  distance  tolerance  and  less  than  five,  the 
score  is  adjusted  by  adding  five  minus  the  distance  to  the  score  [2].  When  all  the  entry 
scores  have  been  generated,  MAK40  returns  to  its  caller  the  completed  S-table  and  the 
number  of  entries  in  the  S-table. 


Distance  = /[(delta  X2  - delta  Xi)^ 

V+ (delta  Y2- delta  Yi)2]  [1] 


Score  = Score  + ( 5 - Distance  ) 


[2] 


To  determine  the  "best  fit"  X and  Y values,  the  routine  WINDOW  is  called.  WINDOW 
generates  values  to  be  placed  in  a table  indexed  by  delta  X and  delta  Y values  from  the  S- 
table.  The  S-table  entries  with  the  same  delta  X and  delta  Y values  are  counted  and  that 
number  is  placed  in  the  proper  delta  X,  delta  Y location.  This  routine  moves  a window  over 
the  delta  X,  delta  Y table  and  generates  three  values  for  each  position  in  that  table.  These 
values  are  for  the  specified  window  size,  a window  size  that  is  two  increments  smaller,  and 
a window  size  that  is  four  increments  smaller.  A starting  window  size  must  always  be  an 
odd  number  and  be  greater  than  or  equal  to  5 units.  This  forces  the  center  of  the  window 
to  be  at  an  integer  coordinate  position,  not  located  between  coordinates.  If  the  specified 
size  is  7x7  then  the  other  two  values  are  for  windows  5x5  and  3x3.  WINDOW  returns  the 
location  of  the  highest  score  for  a large  window.  If  two  or  more  locations  tie,  then  the 
delta  X,  delta  Y position  of  the  one  out  of  this  group  with  the  highest  medium  window 
score  is  returned.  If  a tie  still  exists,  the  highest  score  from  the  smallest  window  has  its 
location  returned.  Finally,  If  a tie  exists  in  ail  three  window  sizes,  the  average  of  delta  X, 
delta  Y locations  that  tie  in  the  smallest  window  are  the  values  that  are  returned. 
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START 
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PROGRAM  FINDTRANS 


O 


I 


READ  IN 
BASE  FINGER 
FILE 


BESTFIT 

BEGIN  LOOP  FOR 
MULTIPLE 
COMPARES 


READ  IN 
FINGER  FILE 
TO  COMPARE 


ASSUME  BEST 
FITATX-0, 

Y-O.e-0 


A-4 


6 


A-5 


1 


BEGIN  S-TABLE 
ENTRY  COMPARI- 
SON LOOP 
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A-6 


A- 7 
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0012 

0013 

0014 
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0016 

0017 

0018 

0019 

0020 

0021 
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0027 

0028 
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0033 

0034 

0035 
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0037 

0038 

0039 

0040 

0041 
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0044 

0045 
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0049 
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0054 

0055 
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0059 

0060 
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0065 
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0070 

0071 

0072 

0073 

0074 

0075 

0076 

0077 
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PROGRAM  FINOTRANS 
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CALCULATES  POSSIBLE  TRANSFORMATION  VALUES 

BY  R.  ALLEN  WILKINSON 

INITALIZE  VARIABLES 

IMPLICIT  INTEGER  (A-Z) 

INTEGER  UNIT. FIN 
INTEGER  NUM. I . J .K. L 

INTEGER  I STAB (500. 6) . I SCNT . CHART (200 . 10) 

INTEGER  BEST(4) .TOTAL(25) . TABLE( 10 . 500 . 6) .TSCNT(10) .GROUP (20) 
INTEGER*2  B(250.3.2) 

CHARACTER  TYPE 

NUM  NUMBER  OF  FINGERS  TO  BE  COMPARED 

ISTAB(K.6)  S TABLE  ENTRIES 
ISCNT  NUMBER  OF  ENTRIES  IN  S TABLE 

CHART(K.L)  ARRAY  FOR  MINUTIAE  MATCHED  CHART 
TOTAL(I)  TOTAL  OF  DISTANCE  CHART  OVER  I FINGERS 
GROUPU)  MINUTIAE  MATCHED  SUMMARY  ARRAY 

DATA  UNIT/7/ 

DATA  TOTAL/0 .0.0. 0.0. 0.0. 0.0. 0.0. 0.0. 0.0. 0.0. 0.0. 0.0. 0.0, 0.0/ 

INCLUDE  ' (WILKINSON. LIB ]COM004.FOR/LI ST’ 

INTEGER*2  A 
INTEGER  IMX 

COMMON  /BLK0O4/A(25O.3,2) . IMX(2) 

A(250.3.2)  MINUTIAE  ARRAY(MINUTIAE  f.DATA  TYPE, PRINT  TYPE) 
DATA  TYPE-  X-1 , Y-2 . THETA-3 
PRINT  TYPE-  SEARCH  PRINT-1. FILE  PRINT-2 
IMX(2)  NUMBER  OF  MINUTIAE(PRINT  TYPE) 

INCLUDE  • [WILKINSON. LIB] COM007. FOR/LI  ST’ 

COMMON  /BLKO07/LN(51 ) . LT(51 ) , LS . LTHETA . MTM2 . I PRMT , UPPER ( 2 ) . 

• ISBZ(2) . ISRZ(2) . I SUN  IT. SCORES (40) . IPRM2 . IJUMP 


1 

c 

INCLUDE 

’ [WILKINSON. LIB ]COM0 12. FOR/LI  ST’ 

1 

LOGICAL 

OIAG 

1 

COMMON 

/BLK012/DIAG(100) 

1 

c 

« • 

FLAG 

SUBROUTINE 

DIAGNOSTIC 

1 

c 

• « 

DIAG( 1 ) 

BESTFIT 

S TABLE  INFO  PRINTED 

1 

c 

« • 

DIAG(2) 

BESTFIT 

S TABLE  DX.DY  VALUES  PLOTTED 

1 

c 

« • 

DIAG(3) 

BESTFIT 

DISTANCE  SUMMARIES  PRINTED 

1 

c 

* • 

DIAG(4) 

BESTFIT 

USING  SCORE  TOTALS  NOT  AVERAGE 

1 C 


1 C 
1 

1 C 
1 C 
1 C 
1 C 
1 C 

1 C 
1 C 
1 

1 C 


1 C 


« • 
« • 
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INCLUDE  ’ [WILKINSON.LIBJCOM013.FOR/LIST’ 


CHARACTER  TMPFILE*40 
INTEGER  DISTAN 

COMMON  /B LK0 1 3/TMPF I LE . D I STAN ( 25 ) 

USED  IN  OPENFBI 

TMPFILE  ORIGINAL  FILE  NAME  UPON  OPENNING 
USED  IN  BESTFIT 

DISTAN(I)  DISTANCE  FROM  CENTER  COUNTS  ON  A FINGER  (I-OISTANCE) 
INCLUDE  ’ [WILKINSON. L I B ] COMDATA . FOR/L I ST ’ 

DATA  ARRAYS  • 

I NT EGER *2  F I NGER „ DUAL  I TY . MI NCNT . XYT 

COMMON  /BLKDATA/FINGER( 10) .QUALITY(10) .MINCNT( 10) , 

• XYT(250.3. 10) 

USED  IN  NEW  FBI  DATA  ORGANIZATION 


1 

c «• 

FINGER(  I ) 

I th 

FINGER'S 

ACTUAL  NUMBER 

1 

C •• 

QUALITY( I ) 

I th 

FINGER’S 

DATA  QUALITY 

1 

C «• 

MINCNT( I ) 

I th 

FINGER’S 

NUMBER  OF  MINUTIAE 

1 

C •• 

XYT(J .K.  I) 

I th 

FINGER’S 

Jth  MINUTIAE  DATA 

1 

C •• 

WHERE  K- 

1 

C •• 

1 FOR  X 

DATA  VALUE  0-255 

1 

C •• 

2 FOR  Y 

DATA  VALUE  0-255 

3 FOR  THETA  DATA  VALUE  0-359 

INCLUDE  ’ [WILKINSON. L I B ] COM 1 0 . FOR/L I ST ’ 

STRUCTURE  /FBIDAT/  a o 

UNION 
MAP 


0084 

1 

CHARACTER*10  PCNUM 

0085 

1 

CHARACTER»40  CLASS 

0086 

1 

CHARACTER*9  CDATE 

0087 

1 

CHARACTER-8  CTIME 

0088 

1 

CHARACTER*3  FILL 

0089 

1 

INTEGER*2  CARD 

0090 

1 

INTEGER*2  OFFSET(10) 

0091 

1 

INTEGER«2  VALUES(7550) 

0092 

1 

END  MAP 

0093 

1 

MAP 

0094 

1 

CHARACTER*15192  BUF 

0095 

1 

END  MAP 

0096 

1 

END  UNION 

0097 

1 

END  STRUCTURE 

0098 

1 

RECORD  /FBIDAT/INPUT 

0099 

0100 

1 C 

1 

COMMON  /BLKIO/INPUT 

0101 

1 C •• 

USED 

IN  NEW  FBI  I/O  ROUTINES  - FB lOPEN . FB I READ , FB 

0102 

1 C •• 

PCNUM  PROCESS  CONTROL  NUMBER 

0103 

1 C •• 

CLASS  CLASSIFICATION 

0104 

1 C •• 

CDATE  FILE  CREATION  DATE 

0105 

1 C •• 

CTIME  FILE  CREATION  TIME 

0106 

1 C •• 

FILL  JUST  USED  TO  FILL  TO  WORD  BOUNDARY 

0107 

1 C •• 

OFFSETfl)  OFFSET  FROM  PCN  TO  FINGER  #I*S  DATA 

0108 

1 C •• 

VALUES(I)  ACTUAL  MINUTIAE  DATA  ARRAY 

0109 

01  1 0 
0111 

1 C 

C 

DIAGI 

1 )-. FALSE. 

01  1 2 

DIAGI 

2)-. FALSE. 

01  13 

DIAGI 

3)-. FALSE. 

01  14 

DIAGI 

[4)-. TRUE. 

01  15 

01  1 6 

C 

C •• 

OPEN 

FILES 

01  1 7 
011  8 

C 

WRITE(6, 1 100) 

01  1 9 

1 100 

FORMATf IX. ’BASE  FINGER’) 

0120 

CALL 

FBIOPEN(UNIT) 

0121 

0122 

C 

C •• 

READ 

RECORD 

0123 

0124 

C 

CALL 

FBIREAD(UNIT) 

0125 

CLOSE(UNIT) 

0126 

CALL 

FINNUM(FIN) 

0127 

IMX( 1 )-MINCNT( FIN)  ISET  NUMBER  OF  MINUTIAE 

0128 

DO  10  1-1  . IMXM  ) 

0129 

DO 

1 5 J-1 . 3 

0130 

A(  I , J . 1 )-XYTf I . J .FIN) 

0131 

B(  I . J . 1 )-XYT( I . J .FIN) 

0132 

1 5 

CONTINUE 

0133 

1 0 

CONTINUE 

0134 

0135 

C 

C •• 

ENTER  LOOP  FOR  ALL  TRANSFORMING  FINGERS 

0136 

0137 

C 

WRITE(6 . 1 1 01 ) 

0138 

1101 

FORMATfIX. ’ENTER  NUMBER  OF  FINGERS  TO  COMPARE  TO 

0139 

READ 

(5. ’ ( 12) • )NUM 

0140 

DO  20  L-1 .NUM 

0141 

0142 

C 

C •• 

OPEN 

FILE 

0143 

0144 

C 

WRITE(6. 1 102) 

0145 

1 1 02 

FORMAT(1X. ’FINGER  TO  TRANSFORM’) 

0146 

CALL  FBIOPEN(UNIT) 

0147 

0148 

C 

C •• 

READ 

RECORD 

0149 

0150 

C 

CALL  FBIREAD(UNIT) 

0151 

CLOSE(UNIT) 

0152 

CALL  FINNUM(FIN) 

0153 

BESTf 1 )-0 

0154 

BESTf 2)-0 

0155 

BEST(3)-0 

0156 

IMX(2)-MINCNT(FIN)  !SET  NUMBER  OF  MINUTIAE 

0157 

DO 

30  J-1 . IMX(2) 

0158 

DO  35  K-1 .3 

0159 

0160 
0161 

35 

A(  J . K . 2)-XYTf J . K . FIN) 

B( J . K . 2)-XYT( J . K . FIN) 

CONTINUE 

0162 

30 

CONTINUE 

0163 

0164 

C 

C •• 

DO 

FINE  ALIGNMENT  CHANGES  TO  FIND  BEST  FIT 

0165 

0166 

C 

LS- 

-30 

A-9 


0167 

0168 

0169 

0170 

0171 

0172 

0173 

0174 

0175 

0176 

0177 

0178 

0179 

0180 

0181 

0182 

0183 

0184 

0185 

0186 

0187 

0188 

0189 

0190 

0191 

0192 

0193 

0194 

0195 

0196 

0197 

0198 

0199 

0200 

0201 

0202 

0203 

0204 

0205 

0206 

0207 

0208 

0209 

0210 

021  1 

0212 

0213 

0214 

0215 

0216 

0217 

0218 

0219 

0220 

0221 

0222 

0223 

0224 

0225 

0226 

0227 

0228 

0229 

0230 

0231 

0232 

0233 

0234 

0235 

0236 

0237 

0238 

0239 

0240 

0241 

0242 

0243 
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910 

900 

C 


« « 
m m 
« « 


LTHETA-1 2 
DIAG(  1 )’ 
DIAG(2 


FALSE. 

FALSE. 


DIAGCsi-. FALSE. 

CALL  BESTFIT(BEST.B) 

LS-15 

CALL  BESTFIT(BEST,B) 

LS*8 

CALL  BESTFIT(BEST.B) 

LS-4 

DIAGf 1 TRUE. 

DIAG(2 TRUE. 

DIAG(3)-. TRUE. 

CALL  BESTFIT(BEST.B) 

DO  21  Z-1 .25 

TOTAL (Z)-T0TAL(Z)+0 1 STAN (Z) 

CONTINUE 

CALL  TRNSFX(A.2. IMX(2) . BEST(3) , BEST(2) .BEST(1 )) 
CALL  MAK40( ISTAB . ISCNT) 

TSCNT( L)-ISCNT 
DO  23  Z-1 . ISCNT 
DO  22  K-1 . 6 

TABLE( L,Z.K)-ISTAB(Z,K) 


NOS 

BEST 

FIT 

FOR 

PRESENT 

LS 

AND 

LTHETA 

VALU 

NDS 

BEST 

FIT 

FOR 

PRESENT 

LS 

AND 

LTHETA 

VALU 

NDS 

BEST 

FIT 

FOR 

PRESENT 

LS 

AND 

LTHETA 

VALU 

NDS 

BEST 

FIT 

FOR 

PRESENT 

LS 

AND 

LTHETA 

VALU 

ADD  CURRENT  DISTANCES  TO  TOTALS 


22 

CONTINUE 

23 

CONTINUE 

20 

CONTINUE 

1000 

WRITE(6. 1000) 

FORMAT( 1 HI  , '•••••••• 

1001 

WRITE(6. 1001 ) 

FORMAT (IX.’  0 1 

2 3 

• . ’ 12  13  14  15 

16  17  18 

WRITE(6, 1002) (TOTAL( 

K) .K-1 .25) 

1 002 

FORMAT( 1 X . 2514) 

C 

C 

INITIALIZE  CHART  VAR 

lABLES 

C 

880 

DO  880  L-1 . NUM 

GROUP( L)-0 

CONTINUE 

895 

DO  890  L-1 . IMX(1 ) 

DO  895  Z-1 . 10 
CHART(L.Z)-0 
CONTINUE 

890 

CONTINUE 

5 6 7 

19  20  21 


8 9 

22  23 


10 

24‘) 


•’) 
1 1 • 


1 004 


930 


920 

C 

C •• 


CREATE  MINUTIAE  MATCHED  ARRAY 

DO  900  L-1 . NUM 

DO  910  Z-1 ,TSCNT(L) 

CHART(TABLE(L.Z,4) . L ) -T AB LE ( L . Z . 5 ) 
CONTINUE 
CONTINUE 

PRINT  MINUTIAE  MATCH  CHART 

COUNT  NUMBER  OF  MATCHES  PER  MINUTIAE 

CREATE  MINUTIAE  MATCHED  SUMMARY 

DO  920  Z-1 . IMX(1 ) 

WRITE(6, 1004)Z, ( CHART (Z. L) . L-1 , 10) 
FORMAT(1X. 14,1014) 

CNT-0 

DO  930  Y-1 , 10 

IF  (CHART(Z. Y) . EQ. 0)  GOTO  930 
CNT-CNT+1 
CONTINUE 

IF  (CNT.EO.0)  GOTO  920 
GROUP (CNT ) -GROUP (CNT)+1 
CONTINUE 

PRINT  MINUTIAE  MATCHED  SUMMARY 


DO  940  Z-1 . NUM 

WRITE(6. 1005) GROUP (Z) .Z 
1005  FORMAT(1X, 13, • MINUTIAE  MATCHED 

940  CONTINUE 
STOP 
END 


13.  ' TIMES’ ) 
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0016 

0017 
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0020 

0021 

0022 
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0025 

0026 

0027 

0028 
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0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 

0058 

0059 

0060 

0061 

0062 

0063 

0064 

0065 

0066 

0067 

0068 

0069 

0070 

0071 

0072 

0073 

0074 

0075 

0076 

0077 

0078 

0079 

0080 

0081 

0082 
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C 

C •• 

C 

C •• 

C 


C 

C •• 
C •• 

C 

C 

C 

C 

C •• 
C •• 
C •• 
C •• 

c •• 
c •• 

c 

c 

c 

1 c 
1 
1 

1 c 
1 

1 c •• 
1 c •• 
1 c •• 
1 c •• 

1 c 

1 c 

1 

1 

1 c 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 c •• 
1 c •• 
1 c •• 
1 c •• 
1 c •• 

1 c 

1 

1 

1 c 
1 

1 c •• 
1 c •• 
1 c •• 
1 c •• 

1 c 
c 
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SUBROUTINE  BESTF I T ( BEST . B ) 


BY  R.  ALLEN  WILKINSON 

INITALIZE  VARIABLES 

IMPLICIT  INTEGER  (A-Z) 

INTEGER*2  B(250,3.2) 

INTEGER  I . J , J1  .K.K1  , VALUE . S , Q . TOT . CNT 
INTEGER  I STAB (500, 6) . I SCNT . SAVE( 4 ) .DIST.NOCNT 
INTEGER  DX,DY.DT.BEST(4) 

REAL  RTEMP 
CHARACTER*8  WORD(3) 

CHARACTER  OUT (-45 : 45 . -45 : 45 ) .OVER *97 


SAVE( I ) 

SAVE  OF 

INITIAL 

BEST 

FIT  VALUES 

0EST(I) 

BEST  FIT 

VALUES 

1-1 

DELTA 

THETA 

1-2 

DELTA 

Y 

1-3 

DELTA 

X 

1-4 

SCORE 

DX 

DELTA  X 

DY 

DELTA  Y 

DT 

DELTA  THETA 

CNT 

COUNTER 

NOCNT 

COUNTER 

SCORES 

NOT 

WORTH  COUNTING 

DIST 

DISTANCE 

BETWEEN 

TWO 

S TABLE  ENTRIES 

DATA  WORD/' ROTATION 

* . 'DELTA 

Y ' 

. 'DELTA  X ’/ 

INCLUDE 

' [WILKINSON 

. LI6]COMe04. 

FOR/LIST’ 

INTEGER*2  A 
INTEGER  IMX 


COMMON  /BLK004/A(250, 3 , 2) . IMX(2) 

A(250,3,2)  MINUTIAE  ARRAY (M I NUT  I AE  # , DATA  TYPE, PRINT  TYPE) 
DATA  TYPE-  X-1 . Y-2 . THETA-3 
PRINT  TYPE-  SEARCH  PRINT-1. FILE  PRINT-2 
IMX(2)  NUMBER  OF  M I NUT  I AE ( PR  I NT  TYPE) 

INCLUDE  ' [WILKINSON. L I B ] COM007 . FOR/ L I ST ' 


COMMON  /BLK007/LN(51 ) . LT(51 ) . LS. LTHETA.MTM2. I PRMT . UPP ER ( 2 ) . 
• ISBZ(2) . ISRZ(2) . I SUN  IT, SCORES (40) . IPRM2. IJUMP 


INCLUDE  ' [WI LKINSON . L 
COMMON  /BLK008/STABLE 
SCUBE(-45 : 45) . 
HEAP(600 .4) , 

IRNG( 1 20) , 
HIGHS(-29:30.3) . 
HWIND.VWIND.AWIND. 
MAXRG1 ,MAXRG2.MAXRNG, 
INCLUDE  • [WILKINSON. L 
LOGICAL  DIAG 
COMMON  /BLK01 2/DIAG( 1 


FLAG 
DIAG( 1 
DIAG(2 
DIAG(3 
DIAG (4 


IB]COM008 . FOR/LIST ’ 

(-45:45.-45:45).  ITABLE  OF  DELTA  Y.X  SCORES 
IPOINTER  TO  HEAP  BY  VALUE  OF  DELTAY 
!UP  TO  300  DELTA  X . THETA . SCORE . FLG/PTR 
I RANGES  FOR  EACH  PLANE 
ICONTIANS  MAX  SCORES  OF  WINDOWS 
{DIMENSIONS  OF  WINDOW 
MAXWIN  {MAXIMUM  WINDOW  VARIABLES 
IB]COM012. FOR/LIST' 

00) 


SUBROUTINE  DIAGNOSTIC 
BESTFIT  S TABLE  INFO  PRINTED 

BESTFIT  S TABLE  DX.DY  VALUES  PLOTTED 

BESTFIT  DISTANCE  SUMMARIES  PRINTED 

BESTFIT  USING  SCORE  TOTALS  NOT  AVERAGE 

INCLUDE  • [WILKINSON. L I B ] COM0 1 3 . FOR/ L I ST  * 


CHARACTER  TMPFILE*40 
INTEGER  DISTAN 


COMMON  /BLK01 3/TMP FI LE. DISTAN (25) 

USED  IN  OPENFBI 

TMPFILE  ORIGINAL  FILE  NAME  UPON  OPENNING 
USED  IN  BESTFIT 

DISTAN(I)  DISTANCE  FROM  CENTER  COUNTS  ON  A FINGER  (I-DISTANCE) 


OVER-'  ’ 

DO  19  Q-1 .25 
DISTAN(Q)-0 
CONTINUE 
DO  20  Q-1  . 4 

SAVE(Q)-BEST(Q)  {SET  SAVE  TO  VALUES  ASSUMED  BEST  FIT 
CONTINUE  {VALUES  FOR  TRANSFORMATION 

BEST(4)-0 


A-11 


0084 

0085 

0086 

0087 

0088 

0089 

0090 

0091 

0092 

0093 

0094 

0095 

0096 

0097 

0098 

0099 

0100 
0101 
0102 

0103 

0104 

0105 

0106 

0107 

0108 
0109 
01  1 0 
0111 
01  1 2 
01  1 3 
01  1 4 
01  15 
01  1 6 
01  1 7 
01  18 
01  19 
0120 
0121 
0122 

0123 

0124 

0125 

0126 

0127 

0128 

0129 

0130 

0131 

0132 

0133 

0134 

0135 

0136 

0137 

0138 

0139 

0140 

0141 

0142 

0143 
01  44 

0145 

0146 

0147 

0148 

0149 

0150 

0151 

0152 

0153 

0154 

0155 

0156 

0157 

0158 

0159 
01  60 
0161 
0162 
0163 
01  64 

0165 

0166 


1 50 


1 80 


40 


50 


70 

52 

51 


55 

53 

60 


181 


1 88 
1 89 


1 95 
1 90 
C 

C * « 


PRINT  •, ’ENTER  BEGINNING  AND  ENDING  VALUES  FOR  ’.WORDCl) 
READ  (5.  • ( 13) ’ )K.K1 
WRITE  (6 . 1 50)WORD( 1 ) . K .K1 

FORMAT  (IHI, ’TABLE  FOR  ’.A8.I4.’  THRU’, 13) 

CNT-1 

WRITE(6. 180)BEST(3) . BEST ( 2 ) . BEST ( 1 ) 

FORMAT  (IX,  ’STARTING  TRANSFORMATION  IS  X- ’ . I 3 . ’ Y- ’ 

• .13.’  THETA-’. 13) 

DO  60  VALUE-K.KI 

DX-0 

DY-0 

DT-0 

TEST-SAVE(1 )+VALUE 

CALL  TRNSFX(A.2. IMX(2) .SAVE(3) . SAVE (2) , TEST ) 

CALL  MAK40( ISTAB , ISCNT) 

S-0 

MAX-1 

NOCNT-0 

DO  40  J-1 . ISCNT 

IF( ISTAB(J .6) . EQ.0)  THEN 
NOCNT-NOCNT+1 
GOTO  40 
END  IF 

DX-DX+ISTAB( J . 1 ) 

DY-DY+ISTAB( J , 2) 

DT-DT4ISTAB( J .3) 

S-S+ISTAB( J . 6) 

CONTINUE 

IF  (ISCNT. EQ.0)  GOTO  70 
IF  ( ISCNT. LT.NOCNT)  GOTO  70 
IF  ( ISCNT. EQ.NOCNT)  GOTO  70 
RTEMP-DX/( ISCNT-NOCNT) 

DX-JNINT(RTEMP) 

RTEMP-OY/( ISCNT-NOCNT) 

DY-JNINT(RTEMP) 

RTEMP-DT/( ISCNT-NOCNT) 

OT-JNINT(RTEMP) 

IF  (ISTAB(J.6) .GT. ISTAB(MAX.6))  MAX-J 
IF  ( . NOT  . DIAGU)  ) then  1DIAG(4)  TRUE  DO  TOTAL 

RTEMP»S/( ISCNT-NOCNT) 

S-JNINT(RTEMP) 

ENDIF 

WRITE  (6,50)DX,DY,DT.S.WORD(1 ) .TEST 

FORMAT  (IX. ’DX: ’ . 13. ’ DY:’.I3,’  DT:'.I3.'  SCORE:’. 15. 

• ’ AT  ’ , A8 , ’ VALUE: ’ . 13) 

IF  (BEST(4) . LT . S)  THEN 

BEST(4)-S 

TOT-TEST 

CNT-1 

ELSE 

IF  (BEST(4) . EQ.S)  THEN 
CNT-CNT+1 
TOT-TOT+TEST 
ENDIF 
ENDIF 
GOTO  51 

WRITE  (6,52)WORD( 1 ) .TEST 

FORMAT  (1X,’ISCNT-0  NO  SCORE  AT  ’.A8.’  VALUE’.  14) 

DO  53  J-1 .250 
DO  55  1-1 .3 

A(J . I .2)-B( J .1,2) 

CONTINUE 
CONTINUE 
CONTINUE 
RTEMP-TOT/CNT 
BEST( 1 )-JNINT(RTEMP) 

WRITE  (6.181)  BEST(3) .BEST(2) .BEST( 1 ) 

FORMAT  (IX. ’BEST  TRANSFORMATION  IS  X-’,I3.’  Y- ’ 

• .13,’  THETA-’. 13) 

CALL  TRNSFX(A.2. I MX ( 2 ) . BEST ( 3 ) . BEST ( 2 ) .BEST(1 )) 

CALL  MAK4@( ISTAB . ISCNT)  IGENERATE  ISTAB 
DO  189  J-1 . IMX(2) 

DO  188  J1-1 .3 

A(J.J1 .2)-B(J . J1 .2) 

CONTINUE 
CONTINUE 
DO  190  J— 45.45 
DO  195  J1— 45.45 
STABLE( J . J 1 )-0 
CONTINUE 
CONTINUE 

A- 1? 

LOAD  PLOT  ARRAY 


0167 

0168 

0169 

0170 

0171 

0172 

0173 

0174 

0175 

0176 

0177 

0178 

0179 

0180 

0181 

0182 

0183 

0184 

0185 

0186 

0187 

0188 

0189 

0190 

0191 

0192 

0193 

0194 

0195 

0196 

0197 

0198 

0199 

0200 

0201 

0202 

0203 

0204 

0205 

0206 

0207 

0208 

0209 

0210 

021  1 

0212 

0213 

0214 

0215 

0216 

0217 

0218 

0219 

0220 

0221 

0222 

0223 

0224 

0225 

0226 

0227 

0228 

0229 

0230 

0231 

0232 

0233 

0234 

0235 

0236 

0237 

0238 

0239 

0240 

0241 

0242 

0243 

0244 

0245 

0246 

0247 

0248 

0249 


C 


201 

202 

200 

C 

C •• 
C •• 
C •• 

C 


240 

250 


252 

251 

C 

C •• 

C 


225 


235 


236 

260 

C 

C •• 

C 

261 

210 


209 

212 


230 


DO  200  J-1  . ISCNT 

IF  ( . NOT . DIAG( 1 ) ) GOTO  202  IPRINT  S TABLE  IF  TRUE 
WRITE  f 6,201  )J  . USTAB(  J ,K)  ,K-1  .6) 

FORMATOx,  • ISCNT:  13,  ’ DX:’.I3,*  0Y;*,I3,’  0T*,I3, 

• ' SM#:'.I3,’  FM#:’.I3,’  SCORE:’, 15) 

STABLER  I STAB (J . 2) , I STAB (J  . 1 ) )- 

• STABLEC ISTABC J . 2) , I STAB (J . 1 ) ) + 1 
CONTINUE 

CONVERT  STABLE  ARRAY  TO  BE  OUTPUTTED 
0-9,  A-Z  , a-z 
0-9, 10-35,36-61 


MAX-0 

DO  250  J— 45,45 
DO  240  J1— 45,45 

IF  (STABLE(J . J1 ) . EQ.0)  THEN 
OUT(J.J1)-’  • 

GOTO  240 
ENDIF 

IF  (STABLE(J . J1 ) . LT. 10)  THEN 

OUT( J . J 1 ) -CHAR ( STAB LE(J , J 1 )+48) 

ELSE 

IF  (STABLEC J . J 1 ) .GT . 35)  THEN 

OUT( J , J 1 ) -CHAR ( STAB LE(J . J 1 )+60) 

ELSE 

OUT(J . J1 )-CHAR(STABLE( J , J1 )+55) 

ENDIF 
ENDIF 
CONTINUE 
CONTINUE 
CALL  WINDOW(5,5) 

BEST(3)-BEST(3)-HWIND 

BEST(2)-BEST(2)-VWIND 

CALL  TRNSFX(A,2, I MX ( 2 ) , B EST ( 3 ) , BEST ( 2 ) , BEST ( 1 ) ) 

CALL  MAK40( ISTAB , ISCNT)  (GENERATE  ISTAB 
DO  251  J-1 , IMX(2) 

DO  252  J1-1 ,3 

A(J,J1,2)-B(J,J1,2) 

CONTINUE 

CONTINUE 

S TABLE  CHART 

IF  ( .NOT.DIAG(2) ) GOTO  261  IPRINT  S TABLE  PLOTTING  IF  TRUE 

WRITE  (6,225) 

FORMAT(1X,’  ....-4 -3 -2 -1  ’ . 

• ’ 0 1 2 3 4*. 

• • ') 

DO  260  J — 45.45 

WRITE  (6.235)J  . (OUT(J.K)  .K— 45,45)  . J 
F0RMAT(1X. I3.91A1 , 13) 

IF  (J.EQ.VWIND)  THEN 

0VER(49+HWIND: 49+HWIND)-’/’ 

WRITE(6.236)OVER 
FORMAT( •+’ . A97) 

ENDIF 
CONTINUE 
WRITE  (6,225) 

SUMMARY  DISTANCE  CHART 

IF  ( . NOT . DIAG(3) ) GOTO  300  IPRINT  DISTANCE  SUMMARY  CHART 

WRITE(6,210) 

FORMAT( 1H1 , ’SUMMARY  CHART  OF  BEST  FIT  DISTANCES’) 

DO  209  J-1 . ISCNT 

D I ST- ( ISTABC J . 1 )-HWIND) **2 
DIST-DIST+C ISTABC J . 2)-VWIND) ••2 
DIST-JNINT(SQRT(FLOAT(DIST)  ) ) 

IF  (DIST.GT.24.0R.DIST. LT.0)  GOTO  209 
DISTAN(DIST41 )-DISTAN(DIST+1 )+1 
CONTINUE 
WRITE  (6,212) 

FORMAT (IX.’  0 1 2 3 4 5 6 7 8 910’ 

• .’  11  12  13  14  15  16  17  18  19  20  21  22  23  24’) 

WRITE  (6.230)(DISTAN(K) .K-1 .25) 

F0RMAT(1X.25I4) 

TEMP-DISTAN( 1 ) +D I ST AN ( 2 ) +D I ST AN ( 3 ) +D I STAN ( 4 ) +D I ST AN ( 5 ) 

PRINT  •, ’TOTAL  0-4: ’.TEMP 
K-TEMP-( IMX( 1 )-TEMP)-( IMX(2)-TEMP) 

RTEMP-( FLOAT (K ))/( FLOAT ( IMX( 1 ) ) ) 

WRITE(6, 183)RTEMP 

A-13  


0250  183 

0251 

0252  184 

0253  300 

0254  182 

0255 

0256 

0257 


FORMAT  (IX.  ’DETECTION  SCORE:  •.F6.3)'®!£i 
WRITE(6. 184) ISCNT. I MX ( 1 ) , I MX ( 2 ) ’ 

FORMAT  (IX. ’ISCNT:  ’.14.’  IMX(1):  ’.14.’  IMX(2):  ’,14) 

WRITE(6. 182)BEST(3) .BEST (2) .BEST(1 ) 

FORMAT  05<.‘BEST  TRANSFORMATION  IS  X-’.I3.’  Y- ’ 

.13.’  THETA-’. 13) 

RETURN 
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SUBROUTINE  FB I OP EN ( UN  I T ) 

OPENS  A FINGERPRINT  FILE 
! ! FOR  VAX  ONLY  1 ! •• 

BY  R.  ALLEN  WILKINSON 

INITALIZE  VARIABLES 


INTEGER  UNIT 

INCLUDE  ' [WILKINSON.LIB]COM013.FOR/LIST' 

CHARACTER  TMPFILE*40 
INTEGER  DISTAN 


COMMON  /BLK013/TMPFILE.O 
USED  IN  OPENFBI 
TMPFILE  ORIGINAL  FILE 
USED  IN  BESTFIT 
DISTAN(I)  DISTANCE  FROM 


STAN(25) 

NAME  UPON  OPENNING 
CENTER  COUNTS  ON  A 


FINGER  (I-DISTANCE) 


OPEN  FILES 

PRINT  •.'ENTER  FILE  IN  SINGLE  QUOTES' 

READ  (5.*)TMPFILE 

OPEN(UNIT , FI LE-TMPF I LE, STATUS- 'UNKNOWN'  . BLOCKS  I ZE-1 5 1 92 . 

• BUFFERCOUNT-1  . RECL-1 5 1 92 . RECORDTYPE- ' VAR  I AB LE ' . 

• ORGANIZATION- 'SEQUENTIAL' . FORM- ' FORMATT ED ' . 

• lOSTAT-IOS , ERR-1000) 

RETURN 

WRITE(6 . 1010) lOS 

FORMAT(1X. '###  ERROR:  '.13.'  ### ' ) 

RETURN 

END 


A-15 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 

001  1 

0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 

0021 

0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 

0058 

0059 

0060 

0061 

0062 

0063 

0064 

0065 

0066 

0067 

0068 

0069 

0070 

0071 

0072 

0073 

0074 

0075 

0076 

0077 

0078 

0079 

0080 

0081 

0082 

A A O ^ 


c 

c •• 
c •• 

c 

c •• 

c 

c •• 

c 


c 

1 c 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 

1 c 

1 

1 c •• 

1 c •• 
1 c •• 
1 c • • 

1 c • • 

1 c * • 
1 c * • 
1 c •• 
1 c 


SUBROUTINE  FB I READ ( UN  I T ) 

READS  A FINGERPRINT  FILE 
! ! FOR  VAX  ONLY  I ! •• 

BY  R.  ALLEN  WILKINSON 


INITALIZE  VARIABLES 


IMPLICIT  INTEGER  (A-2) 

INTEGER  lOS, TEMP, UNIT. CARD. K.PNT 


INCLUDE  ’COMIO. FOR/LIST ' 


STRUCTURE  /FBIDAT/ 

UNION 

MAP 

CHARACTER*10 

PCNUM 

CHARACTER*40 

CLASS 

CHARACTER*9 

CDATE 

CHARACTER*8 

CTIME 

CHARACTER*3 

FILL 

INTEGER*2 

CARD 

INTEGER*2 

OFFSET 

INTEGER*2 

VALUES 

END  MAP 

MAP 

CHARACTER«151 

92  BU 

END  MAP 
END  UNION 
END  STRUCTURE 
RECORD  /FBIDAT/INPUT 


COMMON  /B 
USED  IN  N 
PCNUM 
CLASS 
CDATE 
CTIME 
FILL 
OFFSET! 
VALUES! 


LKIO/INPUT 

EW  FBI  I/O  ROUTINES  - FB I OPEN . FB I READ . FB I WR I TE 
PROCESS  CONTROL  NUMBER 
CLASSIFICATION 
FILE  CREATION  DATE 
FILE  CREATION  TIME 

JUST  USED  TO  FILL  TO  WORD  BOUNDARY 
I)  OFFSET  FROM  PCN  TO  FINGER  #I'S  DATA 
I)  ACTUAL  MINUTIAE  DATA  ARRAY 


1 C 
1 C 
1 

1 C 


INCLUDE  'COMDATA. FOR/LIST ' 

DATA  ARRAYS  • 

I NT EGER *2  F I NGER . QUAL I TY , M I NCNT . XYT 

COMMON  /BLKDATA/FINGERC 10) .QUALITY(10) .MINCNT(10) 
• XYT(250 . 3 . 1 0) 


1 

C 

USED  IN  NEW 

FBI 

DATA  ORGANIZATION 

1 

c 

FINGER! I ) 

I 

th 

FINGER'S  ACTUAL  NUMBER 

1 

c 

QUALITY! I ) 

I 

th 

FINGER'S  DATA  QUALITY 

1 

c 

MINCNT! I ) 

I 

th 

FINGER’S  NUMBER  OF  MINUTIAE 

1 

c 

XYT( J ,K . I ) 

I 

th 

FINGER'S  Jth  MINUTIAE  DATA 

1 

c 

WHERE  K- 

1 

c 

1 FOR  X DATA  VALUE  0-255 

1 

c 

2 FOR  Y DATA  VALUE  0-255 

3 FOR  THETA  DATA  VALUE  0-359 


SELECT  CARD 


WRITE(6 ,•)' ENTER  CARD  TO  BE  READ* 

READ (5 . • )CARD 
C 

C ••  READ  FILES 
C 

10  CONTINUE 

READ(7 . 40 . lOSTAT-IOS , ERR- 1 000 . END- 1 1 00) 

• NBYTES. INPUT. BUF 
40  FORMAT(Q.A) 

IF  (CARD. NE. INPUT. CARD)  GOTO  10 
DO  100  1-1.10 

IF  (INPUT. OFFSET(I) .EQ.0)  THEN 
MINCNT( I )-0 
ELSE 

F I NGER ( I )- INPUT .VALUES! I NPUT . OFFSET ( I )-46+2) 
QUALITY( I )- I NPUT .VALUES ( INPUT .OFFSET! I )-4643) 
M I NCNT ( I )- 1 NPUT .VALUES! I NPUT. OFFSET! I )-46+4) 
PNT- I NPUT. OFFSET! I )-42 
DO  1 50  J-1 .MINCNT! I ) 


0084 

0085 

0086 

0087 

0088 

0089 

0090 

0091 

0092 

0093 

0094 

0095 

0096 

0097 

0098 

0099 


1 50 

1 00 

1000 

1010 

1 100 
1110 

1 1 20 


INPUT . VALUES(PNT+K) 
INPUT .VALUES(PNT+K+1 ) 
INPUT . VALUES (PNT+K+2) 


XYT 
XYT 
XYT 

CONTINUE 
ENDIF 
CONTINUE 
RETURN 

WRITE(6. 1010) lOS 
FORMAT(1X. •###  ERROR 
RETURN 
WRITE(6. 1110) 

FORMAT( 1 X . ' END  OF  FILE  REACHED  1 ! 1 I * ) 
WRITE(6, 1 120)CARD 
FORMAT( IX , 'CARD  # 

RETURN 
END 


13 . ‘ DURING  READ  ### ’ ) 


17 . • MAY  NOT  BE  IN  FILE') 


A-17 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 

001  1 

0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 

0021 

0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 

0058 

0059 

0060 

0061 

0062 

0063 

0064 

0065 

0066 

0067 

0068 

0069 

0070 

0071 

0072 

0073 

0074 

0075 

0076 

0077 

0078 

0079 

0080 

0081 

0082 

0083 


C 

C •• 

C 

C •• 
C •• 
C •• 
C •• 
C •• 
C 


1 C 

1 

1 

1 C 
1 

1 C •• 
1 C •• 
1 C •• 
1C** 

1 C 
C 

1 C 

1 

1 

1 C 
C 

C ** 

C 

C 

C 


SUBROUTINE  MAK40 ( I STAB , I SCNT ) 

PARAMETER ( IPRM1-500) 

ROUTINE  TO  FORM  TABLE  OF  CLOSENESS ...  VALUES 

MAY  22.  1984  1640 

AUGUST  25.  1986 
SEPTEMBER  18.  1986 

TR  CALCULATED  AS  THE  HYPOTENUSE  RATHER  THAN  BY  JUST  LIMITS 
IN  X AND  Y 

IMPLICIT  INTEGER(A-Q.S-Z) 

REAL  TR . TEMP( IPRM1 ) 

DIMENSION  ISTAB( IPRM1 .6) 

INCLUDE  ’COM004. FOR/LIST ’ 

INTEGER*2  A 
INTEGER  IMX 

COMMON  /BLK004/A(250.3.2) . IMX(2) 

A(250.3.2)  MINUTIAE  ARRAY ( M I NUT  I AE  # . DATA  TYPE. PRINT  TYPE) 
DATA  TYPE-  X- 1 . Y-2 . THETA-3 
PRINT  TYPE-  SEARCH  PRINT-1. FILE  PRINT-2 
IMX(2)  NUMBER  OF  MINUTIAE(PRINT  TYPE) 


INCLUDE  ’COM007. FOR/LIST’ 

COMMON  /BLK007/LN(51 ) . LT(51 ) . LS . LTHETA,MTM2. I PRMT . UPPER ( 2 ) . 
ISBZ(2) . ISRZ(2) . I SUN  IT. SCORES ( 40 ) . IPRM2. I JUMP 


SETUP  THE  FIRST  2 STEPS.  WHEN  AN  ENTRY  IS  MADE  IN  THE  ISTAB  ON  ONE 
LINE  REPRESENTING  THE  PAIR  (i.j).  THEN  THE  SAME  ENTRY  WILL  BE  MADE 
EVENTUALLY  IN  THE  SLOT  REPRESENTING  THE  PAIR  (j.i). 


60 


90 
1 00 
C 

C ** 


ISCNT-0 
DO  100  II 
LSL-A( I 1 
LSR-A 
LSB-A 
LST-A 
DO  90 


-1 . IMX( 1 ) 

1 . 1 ) -LS 

1 . 1 ) +LS 

2.1 ) -LS 

2 . 1 ) +LS 
1 . IMX(2) 

IF(A( J1 .2.2) .GE. LSB)THEN 
IF(A( J1 .2.2) . LE. LST)  THEN 
IF(A(J1 .1.2) .GE. LSL)  THEN 
IF(A(J1 . 1 .2) . LE. LSR)THEN 
DIST-(A(J1 . 1 .2)-A(I1 . 1 


I 1 

I I 
1 1 
J1 


NUMBER  OF  SLOTS  USED 
CYCLE  TO  100  FOR  EACH  SEARCH  MIN 
LEFT  EDGE  OF  SRCH  BOX 
RIGHT  EDGE 
BOTTOM  EDGE 
TOP  EDGE 

LOOK  AT  EACH  FILE 
ABOVE  BOTTOM  EDGE 
BELOW  TOP  EDGE 
TO  THE  RIGHT  OF  LEFT 
TO  THE  LEFT  OF  RIGHT 
1 ) ) **2+ 


MINUTIA 


EDGE 

EDGE 


DELTA  THETA 


(A(J1  .2.2)-A(I1  .2. 1 ))**2 
IF(DIST. LE. (LS**2))  THEN 
DELTHE-A(J1  . 3 . 2)-A( I 1 . 3 . 1 ) 

IF( lABS(DELTHE) .GT. LTH ETA) THEN 
IF(DELTHE.GE.0)THEN 
DELTHE-DELTHE-360 
ELSE 

DELTHE-DELTHE+360 
END  IF 
ENDIF 

IF( lABS(DELTHE) . LE. LTHETA)THEN 

ISCNT-ISCNT+1  ! FILE  MIN.  WITHIN  BOX 

IFOSCNT.GT.500)  THEN  ! RAN  OUT  OF  ROOM  IN 
PR  I NT  60  11  J 1 

FORMAT(//1x! '###  MATCH87  ERROR  — EXCEEDED 
• ISTAB—  I 1/J1-'  .214) 

STOP 


ISTAB 


ENDIF 

ISTAB( ISCNT. 1 )-Af J1 . 1 
ISTABflSCNT . 2)-A( J 1 . 2 
I STAB! I SCNT . 3)-DELTHE 
ISTABflSCNT . 4)-I 1 
ISTAB( ISCNT . 5)-J 1 
TEMP( ISCNT)-0 
ENDIF 
ENDIF 
ENDIF 
ENDIF 
ENDIF 
ENDIF 

CONTINUE  - 

CONTINUE  A-18 


2)- 
2 


DELTA 
DELTA 
DELTA 
SEARCH  MIN.  # 
FILE  MINUTIA  # 
TS  SCORE 


X 

Y 

THETA 


COMPARE  THE  LIST  OR  SLOTS  A FORM  INITIAL  SCORES 


0084 

0085 

0086 

0087 

0088 

0089 

0090 

0091 

0092 

0093 

0094 

0095 

0096 

0097 

0098 

0099 

0100 

0101 

0102 

0103 

0104 

0105 

0106 


C 


1 40 
1 50 

1 60 


IF(LS.LT.10)  THEN 
KR-LS+1 
ELSE 
KR-1  0 
ENDIF 
KR-5 

DO  150  12-1  . ISCNT 

DO  140  J2-I2+1 . ISCNT  ! LOOK  O THEM  ALL 

IF(I2.EQ.J2)  GO  TO  140  ! DON’T  COMPARE  TO  SELF 

TR-( ISTAB( 12. 1 )-ISTAB(J2. 1 ) ) • • 2+ ( I STAB ( I 2 . 2 ) - 1 STAB ( J2 . 2 ) ) • • 2 

TR-SQRT(TR) 

IF( (KR-TR) .GT.0)THEN 

TEMPf  12  )-TEMPn2)+KR-TR 
TEMP( J2)-TEMP( J2)+KR-TR 
ENDIF 
CONTINUE 
CONTINUE 

DO  1 60  J2-1  , ISCNT 

ISTAB( J2 .6)-JNINT(TEMP(J2) ) 

CONTINUE 

RETURN 

END 


A-19 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 

001  1 

0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 

0021 

0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 
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0059 
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oooooooooooooooooooooooooooooooooooooo  oo 


SUBROUTINE  TRNSFX(A, IK. IMX . DELTAX . DELTAY , DELTHE) 

««  USED  FOR  ROTATING  AND  TRANSLATING  ABOUT  128  x 128 

«•  THIS  ROUTINE  CONSISTS  OF  A SERIES  OF  TRANSFORMS  WHICH  SUPERIMPOSES 

••  THE  CENTRAL  MINUTIA  OVER  THE  FILE  MINUTIA  AND  REPOSITIONS  AND 
••  RECALCULATES  THE  NEW  COORDINATES  AND  ANGLE  FOR  EACH  IMX(1)  MINUTIA 

••  THE  TRANSFORM  TO  BE  APPLYIED  CONSISTS  OF  4 SEPARATE  ONES,  THE  FIRST 

••  [T1]  TRANSLATES  THE  SEARCH  MINUTIA  OVER  THE  CORRESPONDING  FILE  MINUTI 

••  THE  SECOND  [T2]  TRANSLATES  THE  SUPERIMPOSED  MINUTIAE  TO  THE  ORIGIN;  T 

••  THIRD  [T3]  ROTATES  THE  MINUTIAE  BY  -DELTHE  DEGREES  ( NEG.  SIGN  INDICA 

••  A COUNTERCLOCKWISE  ROTATION);  THE  FOURTH  TRANSFORM  [T4]  RETRANSLATES 
••  THE  SEARCH  MINUTIAE  GROUP  OVER  THE  CORRESPONDING  CENTRAL  FILE  MINUTIA 

••  [X-.YM]  - [X.Y.1]  [T1]  [T2]  [T3]  [T4] 

••  THE  TRANSFORMS  MAY  BE  EXPRESSED  AS: 
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1 

0 

0 
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-SIN 
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RY 

1 

WHERE:  RX0/RY0  ARE  DELTAX  AND  DELTAY  AND 

RX  / RY  ARE  THE  X 4 Y COORDINATES  OF  THE  CENTRAL  FILE  MINUTIA 

••  AFTER  ALL  THE  4 TRANSFORMS  ARE  SIMPLIFIED  (MULTIPLIED)  TOGETHER  THE 
NEW  X AND  Y VALUES  OF  THE  MINUTIA  MAY  BE  EXPRESSED  AS  : 


X ' 

- COS( 

X+X1  ) 

+ SIN( 

Y+Y1  ] 

1 + RX 

Y ' 

- COS( 

Y+Y1  ) 

- SIN( 

X+X1  ] 

1 + RY 

WHERE  XI  - DELTAX-RX  ; AND  Y1  - DELTAY-RY 

TRNSFX(MINUTI AE  DATA  ARRAY . F I NGER , NUMB ER  OF  MINUTIAE 
••  OF  FINGER. CHANGE  IN  X. CHANGE  IN  Y. CHANGE  IN 

••  THETA) 


1 03 


IMPLICIT  INTEGER  (A.Z) 
INTEGER*2  A . X 1 , X2  . Y 1 . Y2 . RX . RY 
INTEGER*2  DELTAX . DELTAY . DELTHE 
REAL  RSIN.RCOS.RAD.RTEMP 
DIMENSION  A(250 .3,10) 
RAD-57.29578 
RSIN-SIN(-DELTHE/RAD) 
RCOS-COS(-OELTHE/RAD) 

RX-1 28 
RY-1 28 

X1-DELTAX-RX 

Y1-DELTAY-RY 

DO  1 03  K-1 . IMX 

X2-A(K . 1 , IK)+X1 

Y2-A(K , 2. IK)+Y1 

RTEMP-RCOS*X2+RSIN*Y2+RX 

A(K . 1 . IK)-JNINT(RTEMP) 

RTEMP«RC0S*Y2-RSIN*X24RY 

A(K . 2 . IK)-JNINT(RTEMP) 

A(K.3.IK)-A(K.3. I K )+DELTHE+360 

A(K  . 3 , IK)-MOD(A(K . 3 . IK)  . 360) 

CONTINUE 

RETURN 

END 


IRADIANS  PER  DEGREE 


ICENTER  X COORDINATE 
ICENTER  Y COORDINATE 
ICALCULATE  X OFFSET  FROM  CENTER 
ICALCULATE  Y OFFSET  FROM  CENTER 


MOVE  BY  X 

OFFSET 

MOVE  BY  Y 

OFFSET 

CALCULATE 

NEW  X 

VALU 

E 
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E 
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NOT 
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0-359 
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0026 

0027 

0028 
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SUBROUTINE  W I NDOW(OFFX . OFFY ) 


• • 
m « 
• • 
• « 


1 c 


1 c 


1 70 
1 60 


1 50 
1 00 


FINDS  HIGHEST  TOTAL  SCORE  IN  A \WINDOW  SIZE 
OFFX  BY  OFFY 

RETURN  COORDINATES  OF  CENTER  IN  VARIABLE 
HWIND  AND  VWIND  (X.Y  RESPECTIVELY) 

IMPLICIT  INTEGER  (A-Z) 

REAL  RTFUP 

DIMENSION  MAX(3) ,SUM(3) 


MAX( I ) 
SUM( I ) 


HIGHEST  SCORE  IN  WINDOId^ 
SCORE  OF  CURRENT  WINDOW 


INCLUDE  • [WILKINSON.LIB]COM007.FOR/LIST* 

COMMON  /BLK007/LN(51 ) . LT(51 ) . LS . LTHETA . MTM2 , I PRMT . UPPER ( 2 ) . 
ISBZ(2) . ISRZ(2) . I SUN  IT, SCORES (40) . IPRM2. IJUMP 

INCLUDE  ’ [WILKINSON . LIB]COM008 . FOR/LIST ’ 

COMMON  /BLK008/STABLE(-45:45. -45:45) , ITABLE  OF  DELTA  Y.X  SCORES 

SCUBE(-45 : 45) . IPOINTER  TO  HEAP  BY  VALUE  OF  DELTAY 

HEAP(600,4).  lUP  TO  300  DELTA  X . TH ETA , SCORE . F LG/PTR 

IRNG(120),  IRANGES  FOR  EACH  PLANE 

HIGHS(-29 : 30 . 3) , ICONTIANS  MAX  SCORES  OF  WINDOWS 

HWIND, VWIND, AWIND.  {DIMENSIONS  OF  WINDOW 

MAXRG1 ,MAXRG2,MAXRNG,MAXWIN  {MAXIMUM  WINDOW  VARIABLES 

MAX( 1 )-0 

MAX(2 )-0 

MAX(3)«0 

DO  100  I — LS,LS  {CENTER  Y 

DO  150  J — LS,LS  {CENTER  X 

SUM( 1 )«0 
SUM(2)-0 
SUM(3)-0 
X-(OFFX-1 )/2 
Y-(OFFY-1 )/2 

DO  1 60  K-I-Y, I+Y  { Y 

DO  170  L-J-X,J+X  {X 

SUM( 1 )-SUM( 1 )+STABLE(K , L) 

IF  (K.GT. (I-Y) .ANO.K. LT. (I+Y) 

SUM(2)» 

. AND . K 


EDGE 

EDGE 


AND . L. LT . 
IP  (K.GT. 
AND. L. LT. 
CONTINUE 
CONTINUE 
IF 
IF 


J + X)) 
I-Y+1 
J + X-1 


LT .MAX( 1 ) 
EQ.SUM(1 ) 


LIMITS 
LIMITS 

{SUM  OF  WINDOW 
AND. L.GT 

-SUM(2)+STABLE(K,  L) 

LT. ( I+Y-1 ) .AND. L.GT 


) SUM(3)-SUM(3)+STABLE(K , L) 


(J-X). 

(J-X+1 ) 


(SUM( 1 ) . 

(MAX( 1 ) . 

IF  (SUM(2) . LT.MAX(2) ) 
IF  (SUM(2) . EQ.MAX(2) ) 
IF  (SUM(3) . LT.MAX(3 
IF  (SUM(3) . EQ.MAX(3 
CNT-CNT+1 
TEMPH-TEMPH+J 
TEMPV-TEMPV+I 


1 50 


GOTO 
THEN 
GOTO 
THEN 
GOTO 
THEN 


SAVE  HIGHEST  SUM 


150 


150 


GOTO 

150 

ENDI  F 

ENDIF 

ENDI  F 

MAX( 1 )-SUM( 

1) 

{ SAVE 

SUM 

MAX(2)-SUM( 

2) 

{SAVE 

SUM 

MAX(3)-SUM( 

3) 

{ SAVE 

SUM 

CNT-1 

{HOW  I 

MANY  AT  THIS  VALUE 

TEMPH-J 

{ SAVE 

Y COORDINATE  OF 

SUM 

TEMPV-I 

{ SAVE 

X COORDINATE  OF 

SUM 

CONTINUE 

CONTINUE 

RTEMP-TEMPH/CNT 

HWIND-JNINT(RTEMP) 

RTEMP-TEMPV/CNT 

VWIND-JNINT(RTEMP) 

RETURN 

END 
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